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Abstract 

We identify the shortcomings of existing ab initio quantum chemistry calculations 
for the hyperfine couplings in the recently characterized azafullerene, C59N. Stan- 
dard gaussian basis sets in the context of all-electron calculations are insufficient 
to resolve the spin density near the cores of the atoms. Using the Projector Aug- 
mented Wave method implemented on top of a standard pseudopotential plane- 
wave density-functional framework, we compute significantly more accurate values 
for the Fermi contact interaction. 
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Electron spin resonance (ESR) is one of the key tools for structural studies of 
defects and radicals. The hyperfine splitting of the resonance line by magnetic 
nuclei is used to identify particular peaks with various sites in the system. To 
make such identifications, a theoretical prediction of the hyperfine splitting 
is necessary. The performance of such predictions has been historically rela- 
tively difficult, because the values of the spin density, which is the required 
observable, are very small and thus hard to resolve. This problem is much 
more pronounced in the cases where the electronic states associated with the 
unpaired electron are p-\ike, which further contributes to the reduction of the 
absolute spin density values near the ionic cores. 

In this paper, we investigate the cause of the especially poor performance 
of the electronic structure calculation for the case of the azafullerene, CsgN, 
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Fig. 1. Structure of azafullerene C59N 



and provide a method that predicts significantly more accurate values for the 
hyperfine splittings. 



The study of doped fullerenes and nanotubes is a very active area, because 
dopants can make fullerenes chemically and electronically active. The potential 
applications range from novel semiconductors to nanomechanical devices and 
even biologically active agents. A general review of heterofullerenes is in [1]. 



In particular, azafullerene, the nitrogen substituted version of Cgo has re- 
ceived considerable attention. The structural and electronic properties of the 
C59N monomer have been studied before [2], and it was found that there is sig- 
nificant derealization of the unpaired electron over the molecular cage. The 
degree of derealization over the C 59 N cage has a significant bearing on the 
electrical properties of the material. However, until recently, C59N was only 
available in dimerized, inactive form. It is now possible to make solid solutions 
of the azafullerene monomer in conventional Ceo i n macroscopic quantities [3], 
and thus direct comparisons with the experiments can be made. The struc- 
ture of the molecule is shown in Figure 1. The molecule has a mirror plane, 
containing the atoms N and C2. 



A calculation of the electronic structure of C59N was published in [3] along 
with the experimental results for the hyperfine splittings. The calculated rel- 
ative values of the hyperfine coupling were showing the same trends as the 
experimental data, but the absolute values were off by more than a factor of 
two. 



At room temperature, due to the rapid reorientation of the molecules in the 
crystal, only the isotropic hyperfine term is observable. Also called the Fermi 
contact interaction, it can be expressed for a particular nucleus by the effective 
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Hamiltonian 

H hS = A bt S I -S e 8(r-R I ), (1) 

where the subscript / refers to the nucleus, and is the isotropic hyper- 
fine coupling constant and S is the spin operator. Because of the Dirac delta 
function, the hyperfine splitting then only depends then on the electronic spin 
density at the position of the nucleus and the nuclear and electronic magnetic 
moments, 

A h{ = -g-7e7/ra(Rj), 

where /i# is the Bohr magneton, 7 e and 7/ are the electronic and nuclear 
magnetic moments, respectively, and h(R) is the electronic spin density at the 
nucleus. The isotropic hyperfine coupling is thus very sensitive to the values of 
the wavefunctions near the atomic cores, which in this case, similarly to other 
conjugated systems, can be very small. There are two sources of error: the 
inaccuracy of the density functional used to describe electronic correlations 
and the incompleteness of the basis set. Due to the aforementioned small 
numerical value of the spin density at the nuclear positions, this latter error 
can be particularly severe. A distinct disadvantage of traditional quantum 
chemistry methods, such as the one used in [3] is the uncontrolled nature 
of the standard basis sets with which the wavefunctions and the density is 
expanded. 

Building on the formal work of Hohenberg, Kohn and Sham [4,5], and the later 
algorithmic work of Car and Parrinello [6] and others [7], it is now practical 
to study large systems with near quantum chemical accuracy using a plane 
wave basis set. Below, we use the traditional density-functional plane-wave 
pseudopotential method and the local spin density approximation (LSDA) 
to relax the molecular geometry, then apply the Projector Augmented Wave 
(PAW) method[8], to calculate all-electron spin densities at the nuclear sites, 
from which the hyperfine couplings may be computed. A key aspect of this 
approach is that PAW is used as a post processing step, thus already existing 
standard plane-wave codes can be used to carry out the computationally in- 
tensive geometrical optimization steps. After that, the hyperfine constants can 
be simply and efficiently computed, even using standard numerical packages. 
The ab initio calculations were carried out with the DFT++ package [9], using 
periodic supercell technique in a 24 Bohr cubic cell; the plane-wave cutoff was 
20 Hartrees. The convergence of the results with respect to the plane-wave 
cutoff were checked by repeating the calculations at 35 Hartrees; the hyperfine 
couplings only changed by less than 10%. We used optimized norm conserving 
pseudopotentials developed by Rappe[10]. 

The PAW method was introduced by Blochl in 1994. It allows the reconstruc- 
tion of all-electron properties from pseudopotential calculations. Originally, 
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PAW was presented as a fully self-consistent all electron method, where the 
reconstructed wavefunctions are used to compute the total energy. Here, we 
just use the method of reconstruction. Accordingly, we define angular momen- 
tum projector functions \p) for every atomic species, which have the property 

where (f>Y( r ) are atomic pseudo wavefunctions and i, j are composite indices of 
the usual angular momentum channels [nlm]. The projectors \p) are localized 
near the nucleus, and they vanish outside the cutoff radius of the pseudopo- 
tential. The reconstruction formula is 

i*™> = i*p*> +£E(i^> - 10) (p R ,^n, (2) 

R i 

where |\l/ ps ) is an extended pseudo wavefunction, ae are the orbitals from an 
all-electron atomic calculation, and R runs through the position of the nu- 
clei. Note that because the atomic pseudo wavefunctions are identical with 
the all-electron atomic wavefunctions outside the cutoff radius, the recon- 
structed wavefunction |\|/ rec ) is identical to the pseudo wavefunction |\I/ ps ) in 
the inter-atomic region. The reconstruction is exact within the frozen core 
approximation if we take a complete set of angular momentum projectors. In 
practice, it is often enough to take just one projector for every pair / and m. 

For our purposes here, we also need the notion of a pseudo operator, which 
arises naturally within the PAW formalism. Using the definition (2), the matrix 
elements of an operator O between reconstructed wavefunctions are given by 

RR',ij 

(3) 

Note that for local and semi-local operators, only on-site terms contribute, 
where R = R'. Thus we can define a pseudo operator O ps , 

o ps = o + £ \p R j) (<^|o|^> - <0 P ^|O|C>) ( PR ,\. 

R,ij 

This pseudo operator, when acting on pseudo wavefunctions, by construction 
will give the same matrix elements as the corresponding all-electron operator 
acting on all-electron wavefunctions. 

To extract the hyperfine splitting from an ah initio calculation, we need the 
all-electron spin density at the nuclear positions, 

h(R) = n 1 (R)-n l (R). 

The simplicity of the delta function operator in (1) made this one of the 
first applications of the above reconstruction ideas. If the angular momentum 
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expansion were complete, the first and third terms of equation 3 would cancel 
exactly, and only the term involving the all-electron atomic states would need 
to be calculated. Also, only the s states have a non-vanishing density at the 
nucleus of an atom, thus (3) reduces to 



where i is the index of the projectors for the s channels, and the weighting 
factors Wi are 



Wi(R) = 




where the subscript "at" refers to densities in isolated atoms. Blochl calcu- 
lated hyperfine constants for a variety doped semiconductor systems [11] using 
equation 4 and obtained satisfactory results by just taking a single projector, 



Equation 5 amounts to a simple rescaling of the pseudo spin density by a 
constant factor. 

The problem with the C59N molecule, as mentioned above, is that the unpaired 
electron is in an orbital which has almost p-like symmetry near the nuclei. 
The unpaired spin density at the nucleus is therefore very small and its value 
becomes very sensitive to how well the charge density is resolved near the 
nucleus. We find that a more robust approach is to retain the projector form 
from equation 3 and use the formula 



Apart from the issue of resolution near the ionic core, the truncation of the 
angular momentum expansion in (3) implies that the projector form (6) is 
expected to underestimate the spin density, while the simple scaling method 
(5) overestimates it by assuming that each s channel has the same scaling 
factor Wi = wo for a given nucleus. Note that this error is in addition to those 
resulting from an incomplete spatial resolution of the wavefunction near the 
core. 

The ionic positions in the C 59 N molecule were relaxed until all forces are less 
than 0.03 eV/A. Figure 2 shows the contours of the charge density in gray, 
and the spin density in black; Figure 3 shows just the spin density. The sig- 
nificant derealization of the unpaired electron is clearly visible, as well as 
the oscillations of the spin density as the distance increases away from the 
nitrogen atom. Table 1 shows our computed hyperfine coupling constants for 



n™ c (R)=Y,n ps (R)w l (R) 



(4) 



n iec (R) = n ps (R)w (R). 



(5) 




(6) 
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Fig. 2. Contour of the charge density (gray) and the spin density (black), of the 
C59N molecule. 




Fig. 3. Contour of the spin density in the C59N molecule. 

selected atoms and the results using the GAUSSIAN all-electron program [12] 
which were published along with the experiment in [3]. It is clear that the 
full projection formula (6) is much more accurate; the error is at most 20% 
compared to experiment. It is interesting that the values obtained by simply 
rescaling the pseudo density are relatively close in agreement with the values 
from the GAUSSIAN calculation. Both traditional approaches overestimate the 
hyperfine coupling by over a factor of two. Note that the experiment only mea- 
sures the absolute value of the hyperfine coupling and that the identification 
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of the experimental values with particular sites is actually inferred from the 
calculated values. Recently, defects in silica were investigated by B16chl[13], 
where PAW was used as an all electron method. There, the conclusion was 
similar to ours: much accuracy can be gained by using the full projection 
method. 

In conclusion, we have computed isotropic hyperfine coupling constants for 
the azafullerene C59N and obtained a much better agreement with recently 
published experimental data than other calculations. We showed that the hy- 
perfine parameters are very sensitive observables in this system because the 
unpaired electron is more or less in a p-like state with relatively small spin 
density at the nuclear sites. Thus, a full projection-based formula of the PAW 
method was needed to obtain reasonable values for the hyperfine coupling, as 
opposed to the currently accepted method of simply rescaling of the pseudo 
density. This underscores the significance of basis set convergence and related 
numerical issues regarding the traditional quantum chemistry basis sets. 
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method 


method 


B3LYP 


N 


0.36 


0.33 


0.62 


0.87 


C2 


1.18 


1.38 


2.57 


2.29 


C3 


0.48 


-0.40 


-0.74 


-0.80 


C4 


0.52 


0.55 


1.02 


0.88 


CIO 


0.25 


0.30 


0.56 


0.44 



Table 1 

Hyperfine coupling constants in mT. The present work is in the second and third 
column. The "PAW projector method" uses equation 6, the "scaling method" uses 
equation 5. Note that the experiment only measures the absolute value, so this is 
what is presented in the bar graph. The GAUSSIAN and experimental values are from 
[3]. 
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